Chaotic filtering of moving atoms in pulsed optical lattices 
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We propose a mechanism for a velocity-selective device, which exploits the fundamental phe- 
nomenon of dynamical localization. It would allow packets of cold atoms travelling through a pulsed 
optical lattice in one direction to pass undisturbed, while dispersing atoms travelling in the opposite 
direction. This mechanism is based on a chaotic diffusion process, with a momentum-dependent 
diffusion coefficient. We analyse this classical effect, and we show that it is enhanced by dynamical 
localization. 
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There is much current interest in the development of 
new techniques to manipulate cold atoms. Recent work 
in atom optics has resulted in new devices termed "atom 
chips" 0,0 where cold atoms may be trapped and guided 
by fields above a solid substrate. One might even in 
future envisage a multi-component chip combining arrays 
of micro-traps and atom wires, forming a sort of atom 
optics version of an integrated circuit. Within such an 
atom chip, techniques to, for example, split, transport 
and otherwise control the traffic of atoms, might play an 
important role. 

Optical lattices in particular are receiving much cur- 
rent attention in atom optics, in a wide range of exper- 
imental studies. For instance, cold atoms in optical lat- 
tices have become a paradigm in the study of quantum 
chaos; experiments on sodium and cesium atoms Q pro- 
vided a textbook demonstration of Dynamical Localiza- 
tion (DL) 0, the quantum suppression of classical 
chaotic diffusion. In these experiments the atoms expe- 
rience a periodically pulsed or driven standing wave of 
light. 

The corresponding classical motion is fully chaotic for 
sufficiently strong driving. This implies that the energy 
of the system grows diffusively with each consecutive 
kick: in the absence of phase-space barriers, which are 
present only in the regular regime, the average energy 
(p 2 ) of the particles is unbounded and this diffusive in- 
crease in energy continues indefinitely. It is characterised 
by a diffusion rate Do, ie (p 2 ) = Dot. However in the 
quantum case this process is suppressed on the time- 
scale of the so-called "break-time" t* ~ D /h 2 . After t* 
the atoms absorb no more energy. 

The series of ground-breaking experiments in [j| was 
followed by other experiments with optical lattices prob- 
ing a wide range of quantum chaos phenomena including 
dynamical tunnelling jfj, |7fl, the effect of quantum loss 
of coherence on dynamical localization § and quantum 
accelerator modes 0. The experiments in Q employed 
BECs in optical lattices. The dynamics of Bose-Einstein 
Condensates in optical lattices has recently become of 
very broad interest [To|. 

Here we show that a straightforward modification of 
the quantum chaos experiments can form the basis for a 
device to control the traffic of cold atoms moving along a 



channel in, say, an atom chip by selecting a specified 
velocity. It was recently proposed [llj that atoms in 
double-well lattices, pulsed with unequal periods, will 
have diffusion rates which are momentum dependent: 
D(jp) = Dq + C{p) where C(p) is a momentum depen- 
dent correction. This can form the basis for a type of 
ratchet dynamics: for a system with zero initial momen- 
tum, (p{t = 0)} = 0, particles with p > 0, say, absorbed 
energy at a different rate from particles with p < 0, lead- 
ing to a non-zero net current. While for that system D(p) 
is locally asymmetric around p = 0, it has no particular 
symmetry over a larger scale and has complex oscilla- 
tions with respect to p, with contributions from several 
terms. Hence it is not suitable for the velocity selector 
we propose here. Below we show that there is a system 
which has a diffusion coefficient of very simple form, with 
a single odd-parity correction term, ie C(p) = —C{—p). 
For experimentally accessible parameters we can then se- 
lect a regime where, for particles moving in one direc- 
tion D + (p) = Do + C(p) ~ which means they ab- 
sorb little or no energy; while particles moving in the 
opposite direction experience an enhanced diffusion rate 
D~{—p) > Dq. While in the classical case, this effect 
is confined to a narrow parameter range, we show here 
that for the equivalent quantum system, this "filtering 
effect" is in fact substantially enhanced by DL relative 
to the classical case, and remains effective over a wider 
parameter range. We investigated DL in this anoma- 
lous diffusion regime and find that it is associated with 
a local "break-time" t*(p) ~ D(p)/h 2 . The ratio of en- 
ergy absorbed by particles moving in opposite directions, 
classically is ~ D + (p)/D~(p) for short times and ~ 1 for 
large t ; in the quantum case, due to dynamical localiza- 
tion, the corresponding ratio is [D + (p) / D~ (p)) 2 ~ est for 
t > t* . If implemented, this would represent an atom op- 
tics application of what to date remains a much studied 
fundamental physics phenomenon. 

In the dynamical localization experiment of 0, 0, @ j 
the dynamics is approximately given by the kicked- 
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particle Hamiltonian : H = ^ — K cos x^2 n 8(t — nT) 
where K is the kick strength. The classical dynamics is 
obtained by iterating the well-known "Standard Map": 
££-l_i = Xi +PiT; pi+\ = pi +K sinajj+i. In the Standard 
Map, we can take T = 1, without loss of generality, but 



for the proposed system, we use a repeating cycle of un- 
equally spaced kicks. For simplicity, we take a length-2 
cycle, with the spacing between kicks alternating between 
Xi and T 2 . The Hamiltonian is now given by : 
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V(x)^2^26(t-nT tot + J2Ti) (1) 



n=Q M=l 



with T to t =Ti + T%. This means that the first kick after 
t = is at t = Ti, the second kick at t = T% +T2, the next 
one at t = T\ + (Ti +T 2 ) and so on. We will consider only 
a small deviation from equally spaced kicks, which can be 
defined with the small parameter b, T\ = l + b,T 2 = 1 — 6. 
The spatial symmetry is broken by a linear term which 
alternates in sign: V(x) = — (Kcosx+Ax(— where j 
is the kick number. This Hamiltonian leads to the map: 

Xi = Xi-! + b) 

Pi = Pi-i + K sin Xi + A 
Xi+i = Xi +Pi(l - b) 
P1+1 = Pi + Ksinxi+i - A 



This classical map was investigated previously |l2j| in the 
regular regime. Here we have investigated the proper- 
ties of this map in the chaotic regime. We have also in- 
vestigated the equivalent quantum behaviour. Since the 
Hamiltonian involves only delta- kicks, we consider the 
usual time evolution operator in a matrix representation 
using a plane- wave basis. We find that while single kicks 
couple different quasimomenta, the combined time evo- 
lution operator for two kicks does not. Hence it is most 
convenient to iterate the two-kick time-evolution opera- 
tor. In the usual plane-wave basis \n) = exp (inx) this 
takes the form of a matrix U(2) with elements: 

(n\U q \l) = e -m+b)(i+ q ) 2 h) 

E e " i(1 ^ (i+<?a)2ft ^-^(f) J i— ( 2 ) 



where q is a quasimomentum, ka — int(q — A) and qa = 
q — A — ka. We solve for the quantum time evolution by 
simply iterating repeatedly ip(2) = Uip{t = 0). In Fig. Q] 
we demonstrate the effect on two quantum wavepackets 
by iterating Eq. J2J over 50 pairs of kicks for two different 
sets of parameters (K = 2.0 and h = 0.25 for (a), K = 3.2 
and h = 1.0 for (b)). We show that the device functions 
as a sort of Maxwell demon: while packets moving to 
the right pass the pulsed lattice relatively unperturbed, 
those moving to the left are strongly dispersed. Below 
we analyse this effect. 

For the Standard Map, at the lowest level of approxi- 
mation, the classical momenta at consecutive kicks are 
uncorrelated and evolve in time as a "random-walk" . 
This means that the average momentum of a large en- 
semble of particles is unchanged, and that the average 
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FIG. 1: Demonstration of the filtering effect: figure compares 
the effect of 50 pairs of pulses on two quantum wavepackets 
moving in opposite directions. Initial wavepackets are shown 
with a thin line, while the final wavepackets are shown with 
a bold line. Parameters : (a) for kick strength K — 2.0, h = 
0.25, (b) for K = 3.2, h = 1.0. The values of A and b are, in 
both cases, tt/2 and 0.01 respectively. The initial wavepackets 
are centred around ±po, with 46po = 7r. The initial widths of 
the wavepackets are similar to the experimental values. The 
figure shows that while the right-moving wavepacket is only 
slightly perturbed, the left-moving one is strongly dispersed. 



energy grows linearly: the average kinetic energy (p 2 )/2 
grows by K 2 /4 at each consecutive kick. In the absence 
of phase-space barriers the energy is unbounded and this 
diffusive increase in energy continues indefinitely. It is 
characterised by a diffusion rate Dq, ie (p 2 ) = Dot so 
for uncorrelated momenta Dq = K 2 /2. However, corre- 
lations between sequences of consecutive kicks give im- 
portant corrections to the diffusion coefficient. The 2- 
kick correlation for example, originates from correlations 
between adjacent kicks but one, (V (xi)V (xi-2)) , and 
gives a correction —K 2 J2(K) to the diffusion coefficient. 
The diffusion coefficient, including the first corrections, 
is: D = K 2 {l/2-J 2 (K)-(J 1 {K)) 2 + ...) [I![l|. These 
corrections have even been measured experimentally with 
cold cesium atoms in pulsed optical lattices . 

For the Hamiltonian in Eq. QJ, the correlations take 
a modified form. The corrections they induce are now 
momentum dependent: for an ensemble of particle with 
initial momentum po, the average energy spread at time 
t is given by 

<(p - pof) = K 2 t (I - \[J X {K{1 + b)) 2 + J X {K{1 b)) 



K 2 $(t)C(2,p ) 



(3) 



Where 



C(2,po) 



Jo{2bK) cos(2pob~ A{1 - 
[J 2 (K(l + b)) 



b)) 

.J 2 (K(l-b))} 



3 



and 



*(*) 



1 - J {2bKy- 2 
1 - J Q (2bK) 2 



Higher order correlations (4- kick, etc.) involve more 
Bessel functions, and induce corrections that may be- 
come significant over part of the parameter range, but 
eq. © has the main features. The first correction in 
equation © is due to 3-kick correlation, and is very sim- 
ilar to the J\(K) 2 correction for the standard map. The 
C(2,p) term is the modified version of the 2-kick corre- 
lation term which for the standard map took the form 
—K 2 J 2 (K). This term has several important properties. 
Firstly, it is p-dependent with the cos(2pb — (1 — b)A) 
term. This means that particles with different momenta 
will absorb energy at different rates. This feature is the 
basis of the filtering effect shown above. Secondly, it 
has a non trivial time-dependence, given by the function 
$(f). Since bK <C 1, one has J (2bK) ~ 1 - (bK) 2 . For 
t <C l/(bK) 2 , has a linear behaviour (~ t/2), and 
the C(2,p) term appears as a correction to the diffusion 
coefficient, as in the standard map. However, for larger 
t, $(£) saturates to the value (bK) 2 /2, and the 2-kick 
correlation does not modify the energy growth anymore. 

All the expressions simplify considerably if we consider 
that b ~ 0.01 is a small deviation from period-one pulses, 
so bA, bK << 1. For times t << l/(bK) 2 , we can write: 



«p-po) 2 } = Kh 



J 2 (K) cos(2pob — A) 



(4) 



Hence we have a local diffusion coefficient of the form 
Dip) ~ D - C(2,p), where D ~ K 2 [l/2 - Ji(K) 2 }. 

Figure [21 shows a comparison of this formula to numer- 
ical results: the average energy spread after 20 kicks is 
plotted as a function of K, for an ensemble of 400,000 
particles with a narrow initial momentum distribution 
around p = 7r/(46), with b = 0.005 and A = tt/2. For 
these values, one has cos(2po& — A) = 1. For clarity, we 
have removed the p-independent contribution Dq. We see 
that the numerical results agree very well with the for- 
mula (@J, showing the J%(K) oscillations. For the larger 
values of K shown on the figure, the effect of the func- 
tion $(£) begins to be noticeable, and formula (0} over- 
estimates the result. The full analytic formula (J3J is also 
shown, with an excellent agreement over the whole range 
oiK. 

Figure|31 (a) and (b), shows the average energy spread 
as a function of the initial momentum po, for two different 
sets of parameters. In each case both the classical and the 
quantum results are shown. For the classical case (dotted 
line on the figure) we obtain, as expected from eq. J3J, a 
cosine behaviour in po, with a period w/b. The quantum 
results (full curve) are similar to the classical one. In 
(a), the number of kicks has been chosen to obtain the 
same maximum energy spreading for the quantum and 
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FIG. 2: Average energy spread as a function of K, for an 
ensemble of 400,000 particles with a narrow initial momentum 
distribution around po = 7r/(46), with b = 0.005 and A = tt/2. 
The Dot term has been removed (see text). The numerical 
results are compared to the simple formula (gj, valid for t <C 
l/(bK) 2 , and to the full analytical formula 0. 
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FIG. 3: (a)-(b): average energy spread as a function of ini- 
tial momentum, in the classical case (dotted curve) and the 
quantum case (full curve). The parameters are: (a) K — 3.2, 
h = 1, 60 kicks for the classical curve, 200 kicks for the quan- 
tum curve; (b) K — 1.7, h = 0.25, 100 kicks for the clas- 
sical curve, 100 kicks for the quantum curve. In both case, 
A — n/2 and b — 0.01. In (a), the number of kicks in the 
classical case has been chosen to obtain the same maximum 
value for the energy spread as in the quantum case. The re- 
sults show clearly that the diffusion coefficient oscillates as a 
function of p, with period n/b, and that the minimum of the 
diffusion coefficient is lowered in the quantum case compared 
to the classical one. 

(c) Average energy spread as a function of kick number, for 
two wavepackets: one initially centred around po = n/(4b) 
(the two lower curves), the other initially centred around 
po — —n /(4b) (the two upper curves). The classical results are 
shown as dotted curves, while the quantum results are shown 
as full curves. One can see that the quantum wavepacket with 
slow energy spread localizes extremely quickly, while the other 
wavepacket takes a much longer time to localize. Parameters: 
K = 3.2, b — 0.01, A = 7r/2 and h = 1 in the quantum case. 
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the classical results; this shows clearly that the quantum 
minimum of the energy spreading is much lower than the 
classical one. 

This quantum effect is due to dynamical localization: 
near a maximum or a minimum of the cos(2p$b — A) 
function, the classical diffusion coefficient is enhanced or 
reduced, respectively. Since the break-time, which gives 
the time needed by the quantum system to localize, is 
proportional to the diffusion coefficient, one sees that 
near a maximum of the diffusion coefficient, the quan- 
tum system takes longer to localize, and thus the quan- 
tum wavepacket absorbs energy for a longer time. On 
the other hand, near a minimum of the diffusion coeffi- 
cient, the quantum system localizes more quickly, and the 
quantum wavepacket absorbs energy for a shorter time. 
More quantitatively, this means that the break time now 
depends on initial momentum: t*(po) ~ D(po)/h 2 . As 
a consequence, if the ratio of the maximum and mini- 
mum classical diffusion coefficient (or energy spread) is 
of the order of D+/D-, then for the quantum system 
we expect a ratio for the energy spread of the order of 
(D+/D_) 2 . In fig. El (a) (K = 3.4) for example, the ratio 
D+/D- is roughly of the order of 3, while the ratio of 
the quantum energy spreads is roughly 9. This enhance- 
ment due to dynamical localization is further illustrated 
on fig. El (c), where the energy of the two wavepackets 



(one initially centred around po = n/(4b) (minimum), 
the other around po = — 7r/(46) (maximum)) is plotted 
as function of the kick number, for the classical and the 
quantum case, for K = 3.4. One sees that the quan- 
tum wavepacket with the slowest energy spread localizes 
very quickly (a few kicks are enough), while the other 
wavepacket takes a much longer time to localize. 

In sum, the filtering device results from a combina- 
tion of two physical processes. Firstly, a momentum- 
dependent diffusion coefficient of a particular simple form 
which enables one to select the velocities of interest by 
adjusting the sign and magnitude of b (in the range 
b ~ 0.05 — 0.01 typically). Secondly dynamical local- 
ization, which enhances the classical effect by localizing 
quickly (slowly) the quantum wavepacket which absorbs 
energy slowly (fast). While this quantum filter works well 
over the range K ~ 1.5 — 3.5 the corresponding classical 
one is effective only for a restricted number of kicks and 
over a narrower range of if ~ 1.5 — 1.8. The parameter 
ranges of K and H (H ~ 0.25 — 1.0) are accessible with 
current experiments. 
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